Supersonic dislocations observed in a plasma crystal 
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Abstract 

Experimental results on the dislocation dynamics in a two-dimensional plasma crystal are pre- 
sented. Edge dislocations were created in pairs in lattice locations where the internal shear stress 
exceeded a threshold and then moved apart in the glide plane at a speed higher than the sound 
speed of shear waves, Ct- The experimental system, a plasma crystal, allowed observation of this 
process at an atomistic (kinetic) level. The early stage of this process is identified as a stacking 
fault. At a later stage, supersonically moving dislocations generated shear-wave Mach cones. 

PACS numbers: 52.27.Lw, 52.27.Gr, 61.72.Ff, 82.70.Dd 



Dislocations are ubiquitous in solids They are essential to understanding such prop- 
erties as plasticity, yield stress, susceptibility to fatigue, fracture, and dislocation-mediated 
melting of 2D solids. Dislocation generation and motion is therefore of interest in materials 
science g the study of ea rt h q ua k es 5 snow avalaIKhes Q, colloidal cry8tals & 2D foams 
6| , and various types of shear cracks |7|, |8| . 

The theory of dislocations uses several approaches. In elastic theory, a dislocation's core 
is treated as a singularity in an otherwise continuous elastic material. Linear elastic theory 
predicts that a gliding edge dislocation cannot overcome the sound speed of shear waves 
Ct, because the energy radiated by a moving dislocation becomes infinite at this speed. 
However, a gliding edge dislocation can move at a particular speed of %/2Ct |9( ; in this case, 
it does not emit any radiation at all and therefore its motion is frictionless. A more detailed 
theory of dislocations should take into account the discreet "atomistic" structure of matter. 
In Ref. jsl, this was taken into account in the dislocation's glide plane, outside, the material 
was considered to be an elastic continuum. Then, the limiting speed for dislocations is the 
Rayleigh wave speed Cr; usually, Cr < Ct- 

Gliding edge dislocations moving at the speed of 1.3CV to I.QCt were observed in atom- 
istic computer simulations [l0|. Dislocations were created with a speed exceeding Ct at a 
strong stress concentration. High shear stress was required to sustain their motion. Dis- 
locations always radiated strongly, even when they traveled at exactly \[2Ct- This was 
attributed to nonlinear effects in the dislocation core. Shear cracks propagating faster than 



Ct were observed in Ref. 



8j . However, to the best of our knowledge, experimental evidence 



that dislocations can move faster than Ct is lacking. 



In regular so 
atomistic level 



ids dislocation dynamics is almost impossible to study experimentally at an 



because of the small distances between the atoms (or molecules), high 



characteristic frequencies, and the lack of experimental techniques of visualizing the motion 
of individual atoms or molecules. 

The most suitable model systems to experimentally study dislocations at an atomistic 
level are plasma crystals [3] . These are suspensions of highly charged micron-sized particles 



in a plasma, which can be strongly coupled 



131 ] . Then their mutual interaction causes them 



to self-organize in structures that can have crystalline or liquid order. The interparticle 
distance can be of the order of 100 /im to 1 mm, characteristic frequencies are of the order 
of 10 — 100 s" 1 , and the speed of sound is of the order of 10 mm/s. In addition, the absence 



2 



of a substantial background medium (as in the case of colloids) allows studies of the full 
dynamics without overdamping. These unique characteristics, plus the direct imaging, make 
it possible to study the complex dynamics of crystalline defects [14j . including dislocation 



nucleation and motion 
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151 ]. all at an atomistic level. 



Our experimental setup was a modified GEC chamber as in Ref. 16j, using similar exper- 
imental parameters. Argon plasma was produced using a capacitively-coupled rf discharge. 
We used 42 W of rf power at 13.56 MHz, with an amplitude of 158 V peak-to-peak. The 
self-bias voltage was —96 V. To ensure that the system was not overdamped, a relatively 
ow pressure of 3 mTorr was used. The neutral-gas damping rate is then accurately modeled 



171 ] by the Epstein expression v = 5N g m g v g (p p rp) 1 , where N g , m g , and v g are the num- 



ber density, mass, and mean thermal speed of gas atoms and p p , r p are the mass density 
and radius of the particles, respectively. With leading coefficient 5 = 1.26 this gave 
v = 0.52 s _1 . 

A monolayer of highly charged microspheres was levitated against gravit y in the sheath 

1 I 

above the lower rf electrode. The particles had a diameter of 8.09 ± 0.18 /im |l7( and a mass 
m = 4.2 x 10~ 13 kg. The monolayer included ~ 4500 particles, had a diameter of ~ 60 mm, 
and rotated slowly in the horizontal plane. 

The interparticle potential for particles arranged in a single plane, like ours, is well 



approximated 



181 ] by the Yukawa potential: U{r) = Q(47re r) x exp(— r/X D ), where Q is the 



particle charge and is the screening length. The monolayer is characterized by a screening 
parameter k = b/Xo, where b is the interparticle spacing. In our experiment, b = 0.69 mm. 
We used the pulse technique of Ref. |19|] , making use of a theoretical wave dispersion relation 



20, |2l|], to measure k = 1.58 ± 0.22 and Q = —15 000 ± 1400e. The average sound speeds 



in the central part of our plasma crystal were measured to be Cl = 24.0 ± 1.7 mm/s and 
Ct = 6.3 ± 0.6 mm/s, for compressional and shear waves respectively. (Based on these 
measurements, the Rayleigh wave speed is Cr ~ 6.0 mm/s = 0.95CV, formally using the 
approach of Ref. |9J]). 

The particles were imaged through the top window by a digital camera. We recorded 
movies of 1024 frames at 29.88 frames per second. The 61.1 x 45.8 mm 2 field of view included 
~ 4200 particles. The particle coordinates x, y and velocities v x , v y were then calculated with 
subpixel resolution for each particle in each frame. 

At our experimental conditions, the particle suspension self-organized in an ordered tri- 
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Fig. 1 

FIG. 1: (Color online) Dislocation nucleation and motion in a 2D crystalline lattice. Maps of 
(a) triangulation of the particle positions, (b) bond-orientational function \iPq\, and (c) vorticity 
(V x v) z , where v is the particle velocity, are shown for four different instants of time: (1) 0.335 s, 
(2) 0.569 s, (3) 0.669 s, and (4) 1.004 s. A pair of dislocations, indicated by arrows, was created 
shortly before (2). 

angular lattice with hexagonal symmetry. This lattice always contained defects, as revealed 
by Delaunay triangulation, Fig.[T](a). A defect is defined as a lattice site where a particle has 
a number of nearest neighbors other than six. Defects are highlighted in Fig. Ufa). Most of 
them form linear chains that constitute domain walls in our 2D crystal. Two nearly parallel 
domain walls are seen in Fig. Ufa). 

An edge dislocation in our 2D crystal consists of an isolated pair of 5- and 7-fold defects. 
During the course of experiment, dislocations in the lattice are continuously generated due 
to the shear introduced by the slow rotation. They move around and annihilate with each 
other or with domain walls. Some of them may remain stationary for some time. In this 
Letter, we show that edge dislocations are created in pairs in lattice locations where the 
internal shear stress exceeds a threshold and then move apart supersonically. 
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To study dislocation nucleation, we evaluate the shear strain in the lattice from the 
bond-orientational function \i/jq\ 22j, shown in Fig. [0(b). For every site in the lattice, 
ipQ = ~ Y^j=i exp(6iGj), where Qj are bond orientation angles for n nearest neighbors. In 
the limit of weak simple shear, the following relation can be used: = 1 — 97 2 , where 7 is 
the shear strain. For weak pure shear, = 1 — 2.25e 2 , where e is elongation, which is the 
measure of pure shear deformation. We derived these relations assuming small deformations 
of an elementary hexagonal cell. Note that is insensitive to uniform compressions, 
rotations, and translations, i.e., the deformations that do not involve any variation of bond 
orientations inside an elementary cell. 

To study dislocation dynamics, we use the maps of vorticity (V x v) z , Fig. [0(c), where 
v is the particle velocity (z denotes the out-of-plane component of vorticity). These reveal 
shear motion 23] and are therefore suitable for visualizing moving dislocations. 

The shear strain in the lattice had a non-uniform distribution, as follows from Fig.QJlb). 
The shear strain was higher (or lower) in two kinds of locations. First, it was high in 
domain walls - the two nearly-parallel bright stripes in Fig. QJb) (or equivalently the chains 
of 5- and 7- fold defects in Fig. [H(a)). Second, a "diffuse" concentration of shear strain 
appeared between the domain walls. We attribute this diffuse strain concentration to the 
differential rotation and shear in the lattice, with two "rigid" domain walls imbedded in it. 

The diffuse shear strain increased with time. When it locally exceeded a certain threshold, 
a pair of edge dislocations was created in that location, Fig. W[2). These dislocations appear 
as bright spots in Fig. [TJ[2b) or as pairs of 5- and 7-fold defects in Fig. [Tf2a), all indicated by 
arrows. Once a pair of dislocations was created, they moved rapidly apart, Figs. []](3)J]J4). 

During the course of our 5.7 min experiment we observed about 30 events of dislocation 
nucleation similar to that shown in Fig. [TJ Two more typical examples are shown in Fig. [21 
The "black" and "white" dislocations have opposite directions of the particle rearrangement 
in their cores. 

We examine the dislocation nucleation and dynamics in more detail. The time evolution 
of the shear strain in the lattice location where a pair of dislocations is generated has several 
stages, as shown in Fig. [3]^a). First, the shear strain builds up gradually in a certain location 
(1^6 1 decreases). Second, when the shear strain in this location exceeds a certain threshold, a 
pair of dislocations is born. Third, the shear stress is rapidly relieved when the dislocations 
separate, and gradually drops to the background value. This cycle then starts over again, 
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Fig. 2 

FIG. 2: More examples of dislocation nucleation similar to that in Fig. Qj4c). In (a) "black" and 
(b) "white" dislocation cores, particles rearrange in the clockwise and counterclockwise directions, 
respectively. 

perhaps in a different location. 

When a pair of dislocations is created, their Burgers vectors are oppositely directed and 
equal in magnitude, within the experimental errors and inhomogeneity of the lattice. There- 
fore, the total Burgers vector is conserved in this process. (The Burgers vector represents 
the magnitude and direction of the crystalline lattice distortion by a dislocation.) 

When the two dislocations separate, they leave a stacking fault between them. It appears 
as a narrow (one lattice constant) band where the lattice structure is distorted from trian- 
gular to nearly square, Fig. [T](3a). This stacking fault has a dynamic nature; the lattice is 
restored to its original state by rapid particle rearrangement that is seen as shear motion in 
Fig. [T](3c). A similar shear motion occurs even ahead of the right-hand-side dislocation in 
Fig.[H(2c). 

Next, we analyze the speed of the dislocations as they move apart. Fig. G^b) shows a 
grey-scale space-time plot of vorticity (V x v) z measured along the dislocations' glide plane. 
Superimposed are the positions of dislocation cores. Note that the dislocation motion is 
not smooth; rather, it is of a stick-and-slip type. The average dislocation speeds are higher 
than CV; however, they decrease as dislocations move. (The slope corresponding to CV is 
shown by a dashed line in Fig. EJ^b)). For comparison, the linear shear disturbance in the 
upper part of Fig. Mjo) (seen as a faint black feature) travels at a speed ~ CV- The average 
speeds of dislocations A and B are respectively Ua = 9.5 mm/s and Ub — 13 mm/s. We 
attribute this difference to the radial variation of the local number density of our crystal. It 
diminished from 2.4 mm" 2 in the crystal center (B moved in this direction) to ~ 1.6 mm -2 



6 



1 



0.5 



(a) 



dislocation A 
• nucleation location 
▼ dislocation B 




0.4 0.6 0.8 1 
time (s) 



Fig. 3 



FIG. 3: (Color online) Pair of dislocations is generated and moves apart, (a) Time evolution of the 
bond-orientational function We use to evaluate the shear strain, which is higher when 

\ipe 1 is lower; for details, see the text, (b) The grey-scale space-time plot shows vorticity (V x v) z 
measured along the dislocations' glide plane. Superimposed are the positions of dislocation cores 
calculated as average positions of the respective 5- and 7- fold lattice defects. Dislocations A and 
B are respectively LHS and RHS in Fig. [TJ 

at its periphery where we still observed dislocations (A moved here). 

Dislocations that move supersonically create a distinct signature, i.e., a shear-wave Mach 
cone, shown in Fig. IU(a). The structure is similar to that observed in Refs. 1(3, 



19 



23j. A 



Mach cone is a V-shaped wake created by a moving supersonic disturbance. Mach cones obey 
the Mach cone angle relation sin \i — C/U, where \x is the cone's opening angle, C is the speed 
of sound, and U is the speed of the supersonic disturbance. From their Mach cone angles, we 
derive the relative speeds of dislocations A and B in Fig. EH U a = 1-7CV and Ub = 1.9CV. 
Since we measured Ua = 9.5 mm/s and Ub = 13 mm/s by defect tracking, we can estimate 
the local sound speeds at the locations of A and B as Ct,a = 5-5 mm/s and Ct.b = 6.7 mm/s, 
respectively. Hence, CV was higher in the crystal center and somewhat declined toward 
its periphery due to the corresponding decline in the crystal's number density. A weak 
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FIG. 4: (a) Shear-wave Mach cone generated by a supersonically moving dislocation B from Fig.[3l 
at t = 1.406 s. (b) Distribution of supersonic dislocation speeds U. 



dependence of Ct on number density (or interparticle spacing) was discussed in Ref. 20]. 
The Mach cones created by a dislocation pair were connected by two parallel fronts of shear 
waves, Fig. [T](4c). This wave front configuration evolved from the particle rearrangement in 
the stacking fault, Fig. H](3c). The Mach cones were composed of shear waves and not of 
compressional waves, because they were excited by dislocations moving faster than Ct , but 



slower than Ct- Note that we did not use any external exciter as in Refs. 19j, l23l . 
there were no particles above or beneath the monolayer as in Ref. 25|. 



24|; also 



Finally, we calculated the distribution of dislocation speeds, Fig. IH(b), using the Mach 
cone angle relation and measuring the opening angles of well-developed Mach cones in maps 
like that in Fig. H](a). From this distribution, we conclude that the average speed of super- 
sonic dislocations in our experiment was (1.95 ± 0.2)CV, where Ct is the sound speed of 
shear waves. 

We thank J. Goree and T. Sheridan for their help in conducting the experiment. 
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